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Abstract 

In this work we compile a few differential equations (ODEs) that arise from the relativistic 
tJ" ' equations in cosmological models that consider the "constants" as scalars functions dependent on 

o : 

' time and they are described as perfect as well as viscous fluids. The general idea of the paper is 

5h , to show how to solve the equations of the models through dimensional techniques (self-similarity) . 

I The results are compared with those obtained by other authors and new solutions are introduced. 
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I. INTRODUCTION. 



The purpose of this work is to study different cosmological models that envisage the 
classical constant as scalar functions dependent on time. To outline the differential equations 
that govern such models and to compare their solution through the traditional method 
(integration of ODEs, this is the way followed by other authors) with the Dimensional 
Method (self-similarity). The latter method will consist in the reduction of the number of 
variables and the obtainment of ODEs easily integrable (see^," and in special' for a review 
of self-similarity in General Relativity) 

In all the studied cases, the equations are very similar, there will only be differences when 
the term p (pressure) is defined, depending on if it is considered a perfect or viscous fluid. 
The modified field equations are: (the "constants" G, c and A are functions on t ). 
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where A represents the cosmological "constant" and we impose that the second Bianchi 
identity is verified: (after raising an index) 

(fii-i*|fl)^ = (^I?+A*|)^ (2) 

as well as the so-called conservation principle (bad so-called^ see ' ) for the energy-momentum 
tensor i.e. the covariant divergence of the stress-energy tensor: 

dtv{T^^) = (3) 

this condition will be considered in some cases since we shall see we that can avoid it. 

1. The line element is defined by: 

ds"" = -c^de + f\t) + ide^ + sin^^d^^) 

[1 — kr^ ^ ' 

we only consider the case /c = 0, here. 



^ we understand that this is only an expresion since in GR the conservation of the energy and momentum 
is only approximate 
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2. The energy-momentum tensor is defined by: 

Tij = {p+p)UiUj -pQij 

where p is the energy density and p represents pressure (to be defined in a generical 
way), always verifying [p] = [p]. 

The paper is organized as follows: In the second section, a model will be studied that 
only envisages the "constants" G and A as variable and whose energy-momentum tensor is 
described by a perfect fiuid. As in all the studied cases, this section will begin outlining 
the equations that govern such model, going next to study its detailed solution through 
the two techniques commented previously. The traditional one, will consist in the math 
integration of the equations (non dimensional method). The other is based on dimensional 
techniques. The latter will be studied in two ways, the first one, that we shall designate 
as the simplest method" (naive), will consist in the dimensional study of the problem 
(model) analyzing its set of governing quantities. This set will lead us to the solution of the 
outlined equations through the application of the Pi theorem. While the other dimensional 
technique that we shall study and that we designate as not so simple method, is a mixture 
of the classic dimensional analysis (obtainment of vr-monomials) with the direct integration 
of the equations. The pi monomials will reduce the number of variables in the differential 
equations and therefore it will bring us to obtain a simpler differential equation (ODEs) 
that will be directly integrated. The efficiency of Dimensional Analysis will enable us to 
provide new solutions even in the cases in which we do not envisage some of the departure 
hypotheses as for example div{T-) = 0. In the third section we shall study a model that 
envisages G as well as A as scalar functions and whose tensor is characterized by a bulk 
viscous fiuid. As in the second paragraph, the classic solution of the problem will be studied 
and it will be compared to the one obtained through the dimensional techniques, providing 
new solutions in this case. In the fourth section, constant G, c and A are envisaged as 
dependent functions on time within a model described by a bulk viscous fiuid. Only in 
this case they are provided dimensional solutions since, at the moment, they are the only 
ones known. In the fifth section a model will be studied that may be understood as a 
generalization of the previous one upon considering in it mechanisms of adiabatic creation 
of matter. In this case, it is observed that the governing equations of the model are reduced 
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to the case above therefore no further commentaries are needed. Finally it will be ended 
with some brief conclusions. 



II. CASE G AND A VARIABLE FOR A PERFECT FLUID. 

Attending to the specifications made in the introduction in this case p (pressure) verifies 
the following state equation: 

p = up oj = const. (4) 

where uj G [0, 1] (i.e. it is a pure number) so that the energy-momentum tensor verifies the 
energy conditions. Under these circumstances the equations are: 



f 
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-p + c^A{t) (5) 



(f'Y 8TTG(t) 2wx 



f2 ^2 

from the expressions (5) and (6) we obtain the following relationship 



Gp' + 3il + Lu)pG^ + pG' + ^ = (7) 

furthermore the following law is taken into account 

dzv{T^') =0 ^ p' + 3{p + p)j = (8) 

Now, we go on to see how this model is solved through two methods, one of them, by 
traditional integration and the other by dimensional approach. 



A. Non Dimensional method. 



Deriving the equation (6) and simplifying with (5) we obtain the following relationship. 
This relationship was used by Lau who reconciled the LNH of Dirac with the GR (see ') 



Gp' + 3(1 + uj)pG^- + pG' + -— = 

J OTT 

regrouping terms and taking into account div{T-) = we obtain: 



(9) 



p' + 3(p + p)y 



'^4 



A'c 
8^ 



G' 

Pg 



(10) 



From the equations (9) and (8) the following equation that relates G with A is obtained. 



G' = (11) 

From all these relationships the following differential equation is obtained which is not 
immediately integrable: combining (6)) with (8) (here we are continuing Kalligas et al's 
work see'') 

deriving with respect to t then we obtain: 



pp"-{p'f = l2n{u + If ^ (13) 



that is to say: 

To integrate this equation Kalligas et al made the following hypothesis on the behavior of 
the function G : G oc t° with a G M (see'') this is an unfounded hypothesis in our opinion 
leading to G = Gt" (where G represents certain constant of proportionality). The obtained 
results are: 

' 127r(c^ + l)2Gt"+2 

etc... (see*'). We believe that this solution, even though correct, under the outlined hy- 
potheses presents certain degree of arbitrariness, since we will know neither the behavior of 
G nor the one of any other quantity in function of the equation of state that we are imposing 
i.e. we should give any value to a and guess" which model is describing i.e. what state 
equation belongs to, furthermore they introduce a new dimensional constant, G, of doubtful 
physical meaning. 



B. Dimensional Method. 



In this section, two dimensional tactics are studied. The first of them, designated as 
"nazwe method" (since it is always the simplest method as well as effective) studies the set of 
governing quantities of the problem and through the application of the Pi theorem a solution 
to the equations except the numerical constant ones is obtained. The second method, the 
one designated as "noi so simplest method", combines the classic dimensional analysis with 
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the direct integration of the differential equations. This method has been developed by Prof. 
M. Castans for ordinary differential equations (see' and'"') 



1. Naive Method. 



With the dimensional method that we have followed we are not forced to impose similar 
condition and obviously the results that we obtain are as good or even better since our 
results depend on the state equation that is imposed (see"^ and^*^). The trick is based in 
integrating the equation p' + 3(p + ]?)y = from which we obtain a dimensional constant 
indispensable for our trifling count. With this constant the set of governing parameters is: 
971 = 9Jl(t, c, A^^). With these quantities in a dimensional base 03 = {L, M, T} the problem 
remains perfectly determined with not further conditions (see" and^"). 



2. Not so simple Method 

In this section we shall study several of the possibilities that may arise. 
a. Considering div{T-) = 0. We note in the equation 



A'c^ G' 



SnG ' G 



(14) 



Al A2 

if the conservation principle for the energy-momentum tensor is taken into account, from the 

part [Al] of the equation we obtain the well-known relationship p = A^^f~^^'^^^\ Regarding 

to the second term of the equation {A2) from it, we can extract a vr-monomial tti = Ac^t^ 

that we may express through the following equality A = where c? G M i.e. it is a pure 

number. This vr-monomial is replaced into the equation in the following way: 

dc^ _ ^ioG /'icN 
Ant^ ~ /3('^+i) ^ ^ 

From this equation we can obtain another vr-monomial, vr2 = fc~^t~^ that we may express 

as / = act where a G M is a numerical constant. With this new relationship we simplify our 

last equation, yielding: 

dc^ AM' 



4vrt3 (act)3(^+i) 

that also reads: 

a = ^ f — t^- (16) 
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whose trivial integration is: 



47r A.. ^ ' 



This result has already been obtained through the simplest technique, except the numerical 
constant 

b. Regardless of dw{T^) = 0. If we do not take into account the condition div{Tij) = 
we shall tackle the problem in the following way: 

,' + 3(.+ l)/j^^+,9^ = (18) 

directly we obtain two pi-monomials: tti = Ac^t^ and 7T2 = fc~^t~^ that we replace into the 
equation: 

that we cannot simplify more. At this time, we note that with these hypothesis we do not 
really get anything important. 

c. Considering the conditions G = f3at" and div{T-) = 0. If we take into account 
the hypothesis imposed by Kalligas et al (see*') on the behavior of G = Gt"' with a G M the 
dimensional treatment would be this: first, we define correctly the (dimensional) constant 
that establishes the proportionality between G and t, G = /^ot" in such a way that jSa has 
the following dimensions: [/3a] = Li^ M'^T'"^'" . As above we shall tackle the equation in 
the same way: 



p' + 3(a; + l)p^ 



A'c^ G" 
P- 



iTiG G 



(20) 



Al A2 

from {Al) we obtain p = A^/"'^'-'^^^^ and from (^42) taking into account tti = Ac^t^ and 
G = Paf^ we obtain the following relationship 

dc^ dc^ 1 



f^^^Y-'t^ (22) 



and therefore / 

as we have seen, through the dimensional method the same results can be obtain but in an 
easier way. 
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d. Considering only the condition G = l3at°'. Under this hypothesis we can also 

outhne the problem without considering the condition div{Tij) = 0. With these suppositions 
the equation to solve is: 

p' + 3(. + l),i. + ,- + g^ = (23) 

if furthermore, tti = Ac^t"^ and 112 = pl3at°'^'^c~'^ are taken into account then: 

{2 + a)gc^ , ^ ( gc^ \ f ( Q(? \ a d(? , ^. 

^ ' +3(^ + l)b^ V+U^ (24) 



where g and d are numerical constants. Simplifying 

(2 + a)q , , aq d 



integrating 

d+4TTg(a+l) 

f = K^t^ f = K^t (25) 
where x = '^^3(^^+i)g^^ ^ constant of proportionality with dimensions [K^^] = LT~^ 



As we see, the recipe is always the same. First, consider the equation to be integrated 
obtaining from it the highest number of pi-monomials, in order to decrease the number of 
variables, in such a way that an easily integrable equation is obtained. 

III. CASE G AND A VARIABLE FOR A VISCOUS FLUID. 

This problem was posed by Arbab (see^^). The basic ingredients of the model are: 
The momentum-energy tensor defined by: 

Tij = (p + p*)uiUj - pgij 

where p is the energy density and p* represents pressure [p] = [p*] . The effect of viscosity is 
seen in: 

p* =p- (26) 

where: p is the thermostatic pressure, H = (f'/f) and ^ is the viscosity coefficient that 
follows the law: 

e = ^7P^ (27) 
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where makes the equation be homogeneous i.e. it is a constant with dimensions and 
where the constant 7 G [0, 1]. And p also verifies the next state equation: 

p = up oj = const. (28) 

where uj G [0, 1] (i.e. it is a pure number) so that the momentum-energy tensor verifies the 
so-called energy conditions. 
The field equations are: 

^ll^^lZ^Jj^f^c^Ut) (29) 



3(I2!.?^|W,,e^AW (30) 
deriving (30) and simplifying with (29) it yields 



p' + 3(^ + i)pH - Qkyn^ + ^ + = (31) 

and at the moment we consider this equation. 

f 

Sv{T,') =0 ^ p' + 3{p + p*)j = (32) 

if we develop the equation (32) we get: 

p' + 3{uj + l)pH - Qkjp'^H^ = (33) 

As in the case above, we shall study the model in two ways, one of them analytic (non 
dimensional) and the other dimensional. 

A. Non Dimensional Method. 

In this section we will mainly follow Singh et al work (see'"'). If we take the equation 
(31) regrouped, we get: 



p' + 3{lo + l)pH - dk^p^/H^ = - 



G StiG 



(34) 



Al 



A2 

if take into account the conservation principle 



p' + 3{uj + l)pH - 9k^p"'H^ = (35) 
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then we solve this equation by solving the equation A2 in (34), in such a way that the 
equation to be solved is now: 



■ G' A'c^ 
PT7 



G ■ ° 

this equation is tried to be solved like this (see^'). A = ^^p- is defined (hypothesis by 
Arbab (see^^) as well as by Singh et al (see^"'), condition that as we shall see, is not necessary 
to impose in the solution through D.A.) and from the equation (30) the following relationship 
is obtained: 8nGp = 3(1 — (3)H'^. Therefore if all the equalities are replaced in the equation 
(36) it yields: 



2 H' p' 



{l-(3)H p 

which is easily integrated. 



(37) 



we get to the equation (35) with all these results 

p' + 3(cj + l)pH - gk^p'fH^ = 
we arrive to the next equation: 

p' + 3Ci{u+l)p'^ -gC^k^p^ =0 (39) 



which has got a particular solution in the case 7 = ci ^ obtaining: 



P(^) = 7^ / «o = (3Ci(^ + 1) - 9k^C!) d'' 
[aot) 



and obtaining from it: 



/(t)=C2t(3(^ + l)-3A;,Ci)(l-7) 



This is the most developed solution reached by Singh et al (see '^) which is slightly different 
from the one by Arbab (see^^). 



B. Dimensional Method. 

We shall explore two dimensional methods in this section. The first one, probably the 
simplest one, has the inconvenience of having to depend on Einstein criterion(see^'^ and 
Barenblatt"), while the second one is more powerful and more elaborated. We shall finish 
showing an equation obtained without having to impose the condition div{T-) = 0. 
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1. Naive Method. 

We summarized here the method addressing to the reference (^^) to know more deeply 
the followed method. In this model the set of quantities and constant dJt remain reduced 
to = d)l{t,c,Ai^,k^) while the dimensional base continues being ?B = {L, M,T,6}. The 
followed dimensional method leads us to the obtainment of two dimensionless vr— monomials 
TTi = ^p{'n'2) where ip represents a unknown function. To arrive to a satisfactory solution i.e. 
to get rid of such function ip, we shall have to take into account the criteria of Einstein or 
Barenblatt. With this criterion we obtain the same results as the already existing ones in 
the literature, but let us say, in a trivial way (for more details see^^). 



2. Not so simple Method 

In this section we shall combine dimensional techniques with standard techniques of ODEs 
integration. With the dimensional method we shall obtain dimensionless monomials, which 
will be replaced in the equations. Thus, the number of variables will be reduced in such way 
that the resulting equation is integrable in a trivial way. We study two cases, the first one 
in which we consider div(T-) = 0, while in the other, as we shall see, such hypothesis is not 
needed. 

a. Considering the condition div{T-) = 0. In this case we shall pay attention into 
the equation: 

p' + 3{iu + l)pH - dk^p-^H^ + P% + ^ = 
taking into account the relationship div{T^) = The following equality is brought up: 



p' + 3{u + l)pH - 9k^p^H^ 



Al 



' G' A'c^ 



G SnG 



(40) 



A2 

The idea is the following: By using D.A. we obtain two tt— monomials, which are replaced 
in the equation, obtaining a huge simplification of it. On the other hand we integrate (Al) 
and {A2), solving the problem completely in this way, now without Barenblatt's criterion. 
Let see. The monomials obtained are: tti = pk^^'''t~ and 1x2 = Ac^t^ i.e. 

^ d 

p = akj ''t-r-^ A 



2+2 



6H 
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where a and d are numerical constants. In a generic way the solution is of the following 
way: p = aky'''ti-^ if we define b = then p = ak^t~^ where a = const. G M then 
p' = —hak^^t~^~^ (paying attention only to the term {Al) of the equation) it yields: 

-hak\t-^-^ + 3(cj + l)ak^^t-^H - 9L, [ak^r^ = (41) 

that simplifying it is reduced to: 

9a(^-^) (/')' - Swt-^ff + bt-^f = (42) 



where w = (u + 1), if define 



t 



D 



(w±(w^ -4ba^-^)^ 
(w ± (w^ - 46a^~^)5 



then, the solution has the following form: 



f = lBt 



D 



(43) 



(44) 



(45) 



where / is a certain numerical constant and 5 is a integration constant with dimensions, 
that can be identified with our result by making B = A^k^. 



Now we shall solve the other term of the equation (the A2). the equation 



SwG 







(36)) can be solved in a trivial way if we follow the next results. If we replace the monomials 
TTi = pk^~'^t~ and 112 = Ac^t^ in such equation the integration of it becomes trivial: 



1 , /r^' 



47rGt3 







G' 



dc^ 



aAirk^ 



b-3 



6-2 



(46) 



where a, d and G M (they are pure numbers). We can also observe that this integral needs 
not be solved since a more careful analysis on the number of vr— monomials that we can 
obtain from the equation, leads us to obtain a solution of the type: 



G — G(^k^, c, t) 



which brings us to: 



b'2±b-2 



G{t) = gk-'cH 



(47) 



This method, as we have seen, is more prepared (strived) and the solution, therefore, finer 
though coincident with the previous one. 
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b. Case in which div(T-) = is not considered. Let see now how we can tackle 
(broach) this problem from the D.A. point of view, without imposing the condition div{T-) = 
0. The base *B as before, is still 03 = {L, M, T} while the fundamental set of quantities and 
constant this time is 9Jt = DJl{t, c, /c^), with these data we can obtain the following monomials 
from the equation 

p' + Siuj + l)pH - 9k,p'H' + + ^ = (48) 

considering that: 

p = akf^t^^ ^ = 4^ (49) 



this two monomials are replaced into the equation, which is quite simplified: 



C dr^ 

-bak'^t-'-^ + 3{u + l)ak'^t-'H - 9k^ {ak'^f-'y + ak'^f-'— - —— = (50) 



simplifying this equation, it yields: 

^/ J 2 ^b-3 

-9a^^-'hH' + 3wH - bt~' + ^ - = (51) 

G 4nak^^ G ^ ' 

that along with the field equations (29) and (30) carry us to the next set of equations. For 
example we see: 

= a^k'^t-' + - 

that we replaced into the equation that we are treating (dealing with), resulting: 

( 3c2 + 3(2 j + G Amk\ G 

that solving it results: 

G = gk-'cH'-' (52) 

where (7 G M represents a numerical constant. We finally observe that as in the previous 
section we could have taken into account the three monomials obtained from the equation 
i.e. 

d r^f''-'^ 
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replace them into the equation 

G' A'c^ 

p' + 3(cj + l)pH - QkyH^ + Pt7 + TT^ = 

G ottG 

and calculate /, arriving at the same solution obtained in the section above i.e. 

/ = IBt^ 

We have proved that it is not necessary to impose the condition div{T-) = since, in this 
case, the same solution is obtained as imposing it. 

IV. CASE G,C AND A VARIABLE AND VISCOUS FLUID. 

This model has been developed exclusively through dimensional techniques, therefore, at 
this time no other solution is known, except the one showed in section 4.2. We can consider 
that this model is a natural generalization of the previous one (see' '). 

Following the Alberch and Magueijo work and therefore all the assumptions considered 
there, we arrive to the requirement that the standard equations still retain their form but 
with G{t),c{t) and A(t) varying (see^'^). 

The field equations that describe our viscous model are (in this section we shall take into 
account the hypothesis made in the above section): 



/" , iff 8vrG(t) , 



p* + c\t)A{t) (53) 



developing equation (2) i.e. calculating (^^T/ + A5l) in this case it yields: 

p' + 3(cu + l)pH - gk^H^ + ^ + p^-A-p = Q (55) 

ottG G c 

besides the so-called conservation principle for the energy-momentum tensor is taken into 
account 

*KT/)^0«,' + 3(, + p.)/:^0 (56) 

Now, if we develop the equation (56) it is obtained: 

p' + 3{uj + l)pH - 9kyH'^ = (57) 
14 



p' + 3(cj + l)pH - 9k^p^H^ 



G 



P- 



4ip 

c 



(58) 



Under these circumstances the equation to be treated is: 

8^ ' G 

We cannot play the same trick as before since in this case the "constant" c i— > c{t) varies 
and therefore the equation are much more comphcated. We shall show here two methods: 
the naive and the elaborated one (not so simple). 



A. Naive Method. 



Since in this case we are considering the "constant" c as variable, the set of govern- 
ing quantities remains reduced to 971 = DJl(t, k^, A^) while the dimensional base is still 
!B = {L,M,T). The obtained solutions in this case are perfectly definite since a single 
vr— monomial is obtained. To see more details about this technique we address to reference^"\ 



B. Not so simple Method. 

In this occasion the equation to solve is (considering the condition div(Tij) = 0, as we 
cannot get rid of it): 



p' + 3{uj + l)pH - 9k^pm^ 



— — 

,41 



A'c^ G' n' ' 



(59) 



A2 



This equation, is very similar to equation (40). We shall follow in this paragraph the tactics 
developed in section 3.2.2. On the one hand, we already know how to tackle the left side 
of equation (59) i.e. the Al side. If we take into account these results and imposing the 
following condition A = (that might be unfounded) we can simplify this equation: 



C2(t)t2 

-br^ + 3wH - 9a^^-^kH^ + 



G' 2dc [c't + c] t 



b-3 



G 



- 4- = 

Snak^. G c 



with the usual notation that we are following. From one of the field equations expression 
for G has been obtained. 



G 



Mak^^t-^ 



r, d 



also / is a previous result from last section: 

/ = IBt^ 
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We should not forget that (Al) from equation (59) is the same as the one studied in section 
3.2.2. (40) and p makes no difference. With this expressions we continue our simphfication. 
S.ce . ^ then 

G=<5^1^<- (60) 



We observe again that c = c{t). 



^, ^ 2i3D' - d)cc' {3D'-dy ^,^, 
— -2- 

with this results we go on to simplify the term {A2) in equation (59) obtaining 

G' 2dc[c't + c\t^-^ _ c' _ 



c t 



-2d 



_{3D^-d) 



c t I c 



d (h ^ bd 



2 t 



c = K^t^ (61) 

'h ^ bd 
2~ 6D'\ 

dimensions [K^] = LT~^~^. From this result the expression of G is obtained: 



where x= \ 7; ~ ^ ^ 7, — 2 ) ' proportionality constant that relates c with t with 



With this solution, it is observed for the case 7 = | what corresponds to 6 = 2 

- = = const. (63) 

i.e. the covariance principle is verified. In particular if 6 = 2 then x 
While the cosmological "constant" yields = ^ 



3D^^ 



d d . . 



In this way the remaining quantities are calculated. 
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V. CASE OF G, C AND A VARIABLE WITH ADIABATIC MATTER CREATION. 



The momentum-energy tensor is defined by: 

Tij = {p + P* + Pc)uiUj - {p* + Pc)9ij (65) 

where p is the energy density and p* represents pressure [p] = [p*] . The effect of the viscosity 
is given by: 

p* =p- 3^H 

where: p is the thermostatic pressure, H = (f'/f) and C, is the viscosity coefficient that 
follows the law: 

e = k.^p"' (66) 

The field equations are as it follows: 



ifr 8nG{t) 
/2 3c2(t) 



p + c\t)A{t) (68) 



n' + 3nH - V = (69) 

where n measures the particles number density, ip is the function that measures the matter 
creation, H = f'/f represents the Hubble parameter (/ is the scale factor that appears in 
the metric), p is the thermostatic pressure, p is energy density and Pc is the pressure that 
generates the matter creation. 

The creation pressure pc depends on the function tp. For adiabatic matter creation this 
pressure takes the following form: 

' P + P 



Pc= ~ 



3nH 

The state equation that we next use is the well-known expression 



(70) 



p = ujp (71) 

where u = const, u G [0, 1] physically realistic equations, making in this way that the 
energy-momentum tensor Tij verifies the energy conditions. 
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We need to know the exact form of the function i/j , which is determined from a more 
fundamental theory that involves quantum processes. We assume that this function follows 
the law: 

iP = 3(3nH (72) 

where /3 is a dimensionless constant (if /3 = then there is no matter creation since ip = 0) 
presumably given by models of particles physics of matter creation. 
The generalized principle of conservation brings us to: 

p + 3{uj + 1)(1 - l3)pH - gk^p^H"^ = - 

The conservation principle for the momentum-energy tensor is expressed through the fol- 
lowing law: 

p' + 3{p + p + p,-3^H)H = 

p' + 3{uj + l){l- /3)pH -9k^p^'H^ = (74) 

This equation has been solved in the above section since the only difference between this 
one and the one exposed there is the term (1 — /5) and the solution will be very similar. The 
dimensional solution can be found in the reference (^'). 

VI. CONCLUSIONS. 

The purpose of this work has been to show how to solve this type of models through 
dimensional techniques. We have been able to prove that this technic is applied in a com- 
paratively easy way enabling us to obtain almost trivial solutions and allowing to avoid 
hypotheses needed when using other methods. 

Acknowledgement 1 We wish to express our gratefulness to Prof. Manuel Castans for 
his help and critical commentaries that they have helped the successful attainment of this 
work and to Javier Aceves for his collaboration in the translation into English. 



A'c^ G' , d 



(73) 



BirkhofF, G. Hydrodynamics. -A Study in Logic, Facts and Similaritude (1950) 1st Ed. Prince- 
ton. Sedov, L.I. Similarity and Dimensional Methods in Mechanics. (1959) Infosearch Ltd. 

18 



London. Barenblatt, G. I. Similarity, self-similarity and intermediate asymptotics. (1979) 
Consultants. Bureau New York. Bluman, G. and Cole, J. D. Similarity Methods for Differ- 
entials Equations. (1974) Appl. Math. Sci. N13 Springer- Verlang New York. Bluman, G. and 
Kumei, S. Symmetries and Differential Equations. Springer- Verlang (1989). 

^ Barenblatt, G.I. Scaling, self-similarity and intermediate asymptotics. Cambridge texts in 
applied mathematics N 14 1996 CUP. 

3 Carr, B.J. and Coley, A. A. gr-qc/9806048 

^ Landau&Lifshitz. The classical theory of fields (1971) sec. 101 pp 304 

5 Lau, K-Y. Austr. J. Phys. 38, 547-53, (1985) 

^ Kalligas D, et al. Gen. Rel. Grav. 24, 351-57, (1992) 

^ Castans, M. Davila, P.. XLI Jornadas GIAD Julio 1999 Lanzarote Esparia. 
* Castans, M., Davila, P.. XLII Jornadas GIAD Diciembre 1999 Oviedo Espaha. 
9 Belinchon, J. A. physics/9811017 
° Belinchon, J. A., physics/0002046 

1 Arbab I. Arbab. Gen. Rel. Grav. 29, 61-74, (1997). 

2 Singh, T, et al. Gen. Rel. Grav. 30, 573, (1998) 
^ Einstein, A. Ann. Phys. 35, 679-694, (1911) 

^ Belinchon, J. A., physics/0002022 
^ Belinchon, J. A., physics/9812020 

^ Albrech, A. and Magueijo, J. Phys. Rev. D59, 043516 (1999), Barrow, J.D. Phys. Rev. 
D59, 043515 (1999) 

Belinchon, J. A., gr-qc/9812049 aceppted by Gen. Rel. Grav. 32 N8 (2000) and gr-qc/9907003 
aceppted by Inter. Jour. Theor. Phys. 



19 



